Temperature Dependence of Hall Response in Doped Antiferromagnets 
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Using finite-temperature Lanczos method tlie frequency-dependent Hall response is calculated numerically 
for the t-J model on the square lattice and on ladders. At low doping, both the high-frequency and the d.c. 
Hall coefficient _Rh follow qualitatively similar behavior at higher temperatures: being hole-like for T > Ts! ~ 
!1.5J and weakly electron-like for T < Tg. Consistent with experiments on cuprates, i?H changes, in contrast 
to i?H, again to the hole-like sign below the pseudogap temperature T*, revealing a strong temperature variation 
for T ^ 0. 
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The anomalous behavior of the Hall constant i?H in the nor- 
mal state of cuprates [|l]] remains the challenge for theoreti- 
cians for over a decade. Two aspects, possibly interrelated, 
are evident and should be understood: a) the d.c. i?^ at low 
temperatures T ^ is clearly doping dependent. In the 
prototype material La2-2;Srj;Cu04 (LSCO) it changes from 
positive oc at low doping x < a;*«0.3, consistent 
with the picture of hole-doped (Mott-Hubbard) insulator, to 
the electron-like R'^ < at x > x* in agreement with the 
usual band picture, b) R'^ is also strongly temperature depen- 
dent, both at low doping and optimum doping. At optimum 
doping, the attention has been devoted to the anomalous vari- 
ation of the Hall angle 0^ oc T'^ in YBaaCugOT [§]. On the 
other hand, at low hole concentration Ch < 0.15, Rh{T) in 
LSCO has been shown to follow an universally behaved [||] 
decrease with T in which R'^{T-^0) and the characteristic 



temperature T* of vanishing i?|jT*) ^ both scale with Ch- 
In underdoped cuprates, the same T*{ch) has been in fact as- 
sociated with the (large) pseudogap crossover scale in uniform 
susceptibility Xo(r), in-plane resistivity p{T), specific heat 
c^(T), and some other quantities [^. 

A number of theoretical investigations have addressed the 
first question, i.e. the doping dependence of Rn in models 
of strongly correlated electrons, in particular within the t-J 
model and the Hubbard model on a planar lattice. The advan- 
tage is that one can study the dynamical Hall response and the 
d.c. Hall constant as a ground state (T — 0) property, in partic- 
ular in systems with finite transverse dimension |^, H and in 
the ladder geometry [0, ^. It has been also shown that within 
the t-J model the change from a hole-like to an electron-like 
Hall response can be qualitatively reproduced by studying the 
high-frequency i?y — R^[uj^oo) ijlTI], analytically tractable 
at T ^ oo. Recently, a connection of the reactive i?j^(T=0) 
to the charge stiffness has also been found [^. 

The anomalous temperature dependence of Ru{T), being 
the main subject of this work, has been much less clarified 
in the literature. The Hall mobility HuiT) of a single charge 
carrier in the Mott-Hubbard insulator has been first evalu- 
ated within the generalized retraceable path approximation 
JlO|]. The high-frequency R^(T) has been calculated using 



the high-T expansion [pl]|. At low doping, c/i < 0.15, it has 
been observed that on decreasing temperature is also de- 
creasing instead of approaching presumed (larger) semiclassi- 
cal and experimentally observed d.c. result R^ — l/cheo « 
4i?Jj(r = oo). Related are the conclusions of the quantum 
Monte-Carlo study of the planar Hubbard model [p^, where 
close to the half-filling electron-like R^ < has been found at 
low T. The same has been claimed generally for i?H(^) even 
for low u; [|l2|. Quite controversial are also results for R^{T) 
on ladders In regard to that, we should also mention the 
questionable relation of the off-diagonal axy to the orbital sus- 
ceptibility Xd l§ |l3|], potentially useful as an alternative route 
to the understanding of R^{T) [|l4[. 

In the following we present numerical results for the dy- 
namical Ru{llj), as obtained within the low doping regime of 
the t-J model using the finite-temperature Lanczos method 
(FTLM) [jl^ |l6||. The aim of this letter is to approach the 
low-cj and low-T limit as much as possible and to investigate 
the relation between R^{T) and R^{T). We find these two 
quantities essentially different for T < T*, establishing the 
pseudogap scale T* < J both in the ladder and planar sys- 
tems. 

We study the t-J model in an external homogeneous mag- 
netic field B — curl A, 



H.c. 



F(A)=-t^(e'^^^4d 

+ J ^{Si-Sj - jn.Tij) 



(1) 



where the (inhomogeneous) vector potential enters the phases 
Oij — eA(ri) - Yij. The hopping is only between the nearest 



neighbors {ij). Projected fermionic operators c^^, cj^ do not 
allow for the double occupancy of sites. 

In order to calculate the dynamical Hall coefficient 



Rn{uj) = 



dB 



(Txy{i^) 



B^o Bcrxx{i-^)ayy{uj) 



(2) 



the conductivity tensor is evaluated within the linear response 
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theory, 



le 



Nlu+ 







(3) 



where in the presence of S 7^ the particle current j and the 
stress tensor r operators are given by 



j = *I]ry(ie'«-45^.,+H.c.), 



(4) 



On a square lattice with N sites and periodic boundary con- 
ditions (b.c.) one cannot apply arbitrary magnetic field B 
since only quantized B = B^ — m^oa^/N can be made 
compatible with the periodic b.c. [p^. Therefore the smallest 
but finite B = Bi is used in calculations. The square lattices 
used are in general Euclidean (tilted) N — P + r? , in particu- 
lar we investigate systems N = 10, 16, 18. On the other hand, 
the ladder geometry of N = L x M sites with the periodic 
b.c. in the L direction and open b.c. in the perpendicular M 
direction allows for any finite B ^ Q, the fact already used in 
several T ~ calculations |^, ^, |^. The advantage of lad- 
der systems is also the existence of the reference ground-state 
results i?2(T=0) which seem to be better understood 
Furthermore, at low doping they reproduce the simple semi- 
classical behavior _R2(T=0) ^ — l/chCQ. 

Dynamical components (TQ/3(aj) are evaluated using the 
FTLM [[T^, employed so far for various dynamic and static 
quantities within the t-J model [p^, among them also the 
B = optical conductivity = aaai^^) on a square lat- 
tice. Comparing to the diagonal (Jaa, the evaluation of the 
off-diagonal axy{(^) is more demanding for several reasons: 
a) the introduction of > in the model ([l|) breaks the 
translational invariance and prevents the reduction of the ba- 
sis states in the Lanczos procedure, hence available finite-size 
systems are somewhat smaller, b) we expect (Txij{i^) oc B 
while (Txy{B = 0,0;) does not vanish identically within the 
FTLM; consequently larger sampling over initial wavefunc- 
tions [|l5[ |l^ are needed to reduce the statistical error, c) on 
a finite square lattice the reference result i?[^(r = 0) is not 
meaningful for B„i > 0, while in ladder systems it is quite 
sensitive to the introduction of an additional flux . Never- 
theless, in general, restrictions for the validity of the FTLM 
results are similar to other quantities. Through the thermo- 
dynamic partition function Z(rfs) = Z*, we can define the 
marginal finite-size Tfs below which too few levels contribute 
to the average and results loose the thermodynamic validity 
JT^]. In the following, we analyze results for J = QAt at low 
hole doping Ch — Nh/N {Nh = 1, 2). In this regime we can 
estimate Tfjt ^ 0.15 - 0.2 < 0.5J/t. 

Let us first present results for the dynamical Rn{uj). In 
Fig. 1 we show the normalized real part th = eoCh Re Rh 
for systems with a single hole Nh = 1. In the evaluation of 
Rui^^) from Eq. we insert complex cit B — and 




Figure 1: Dynamical Hall response rn{i^) ~ eoCh Re _Rh(i^) for 
different temperatures T/t and various systems with a single hole 
Nh = 1: a) 2-leg ladder with L = 11, b) 4-leg ladder with L = 4, 
and c) square lattice with A'^ = 18 sites. 



the most sensitive quantity remains axy {(^) calculated at i? = 
Bi on a square lattice and B ^ 0.3i?i on ladders. In the 
presentation of results an additional frequency smoothening 
S = 0.2t is used. The normalization of Ru is chosen such 
that at low doping th = 1 would show up in the case of the 
semiclassical result. 

In Fig. 1 several common features of Ru in the ladder ge- 
ometry and in the 2D systems are recognized: 

a) ru{(^) is quite smoothly varying function of uj, at least in 
contrast to strongly a;-dependent ReCT(a;) on a 2D system, 
which is found [ p^ to decay with an anomalous relaxation 
rate 1/t{uj) (xuj + £,T. 

b) At high temperatures T > t we get a hole-like rn > for 
all systems. In this regime rii{uj) is very smooth, in particular 
for the M = 4 ladder and the 2D lattice. 

c) For low temperatures T < t, ru{uj) is less smooth and the 
dependence is more pronounced for the 2-leg ladder On the 
other hand, M ~ i ladder clearly approaches the behavior of 
the 2D system, whereby both of the latter show quite a mod- 
est variation of rH(w). In all systems the resonances (and the 
variation) visible in ruioj) at high oj > t reflect the predomi- 
nantly local physics of the hole motion and are thus not related 
to a current relaxation rate deduced from a{uj). 

Results for rH(w) are the basis for the calculation of high- 
frequency — rii{uj — oo) as well as the d.c. limit = 
rii{o-!^Q)- The latter is more sensitive since in a finite system 
(even at T > 0) CTq/j (cj-^0) can be singular due to the coherent 
charge transport in a system with periodic b.c. The coherent 
transport shows up in a finite (but small) charge stiffness JT^], 
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Figure 2: D.c. Hall constant rg and the infinite-frequency vs. T/t Yigm& 4: and vs. T/t for different square and ladder lattices 
for various ladders L x M with a single hole Nh = 1. ^ith jv sites and Nh — 2. 
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Figure 3: and vs. T/t for different square lattices with A'^ sites 
and a single hole Nh = 1. 

which should be omitted in the evaluation of Eq. (^. In any 
case, one should take into account proper uj ^ behavior 
of dissipative systems at T > which is different in ladders 
and in 2D lattices, respectively: a) On a ladder we get in the 
leading order of w a normal conductance along the x- 
direction, i.e. (Txx{oj ^ 0) ^ ao, but a finite polarizability 
along the y-direction, ayy{uj^O) oc (^Xyy Hence, we expect 
axy oc u! and finite r^. b) For a macroscopic isotropic 2D 
system we get aaa (i^^O) (Tq and we expect as well a-^y 
(T^ , leading to finite r^. 



In Fig. 2 we present results for r^{T) and r'^{T) for the 
ladder systems with Nh = 1. Results are shown for 2-leg 
ladders with various lengths L — 9, 10, 11 and for M — 4 
ladder with L ~ A. Since and rj!j are properly scaled, 
for given A I curves are expected to approach a well defined 
macroscopic limit at L ^ cxj. In fact, are nearly indepen- 
dent of L (as well as of M) down to T ~ Tfj. A crossover 
at Ts ~ 0.6t from a hole-like rjlj > into a electron-like 
< can be explicitly observed, results are more size 
(L) dependent, nevertheless they reveal a crossover nearly at 
the same T ^ Tg. In contrast to which remains nega- 
tive for the whole regime T < Ts, changes sign again at 
T=T* ^ Q.2t. Although our data for T* are more scattered 
the crossover into the hole-like r^{T<T*) > is expected. 
Namely, from the ground state calculations in same systems 
||] we know that rO(T = 0) - 1.5 and r° (T = 0) - 1.2 for 
M = 2 and M = 4 ladders, respectively. Therefore, it is not 
surprising that the observed dependence r^{T <T*) is very 
steep. 

Corresponding results for the planar lattice in Fig. 3 are 
both qualitatively and quantitatively similar. Note, that at low 
doping the limiting value r'^{T—>oo) = 1/4 agrees with the 
analytical result [pi]], while obtained r'^{T^oo) ~ 0.3 is also 
quite close. Again, the crossover into an electron-like regime 
appears at Tg ^ 0.6t. For larger sizes TV > 16 the lower 
crossover T* ~ 0.2t is visible as well. In finite 2D systems a 
reference numerical result at T = does not exist, however, 
the analytical theory [|l7| indicates that in a macroscopic limit 
with a single hole (iV/, = 1) in an ordered antifeiTomagnet 
one should get = 1. 

In numerically available systems, Nh = 2 represents al- 
ready a substantial doping. Therefore, results for and 
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shown in Fig. 3 should be interpreted in relation with the 
corresponding finite doping Ch- Main message of Fig. 3 is 
that upper crossover Tg, still nearly the same in both r^{T) 
and r^{T), shifts down quite systematically with increasing 
Cfi, i.e. with decreasing size N at given Nh- At least in 
ladder systems at Ch < 0.3, we still find r'^{T = 0) > 
in the ground state ||^], therefore also the lower crossover 
T* < Ts is expected. However, we cannot detect such a 
crossover in r^{T) down to Tf, ^ 0.15t, not surprisingly 
since also the experimental value, e.g. in LSCO at Ch > 0.1, 
is T* < 600K - 0.15i (assuming t - 0.4 eV). 

Let us finally comment on the relation of the d.c. a''^y to the 
orbital susceptibility Xd in a macroscopic 2D system. Namely, 
a°y = eBdxd/dn = eB{dxd/dch){dch/dn), (where /i de- 
notes chemical potential) was derived using seemingly quite 
general thermodynamic relations [||, |l3|], but at the same time 
put under question [|||. Since the d.c. cr°|^(T) > is quite 
a smooth function the above relation seems to yield also a 
qualitative connection between Xd{T) and R^{T). The sit- 
uation should be particularly simple at low doping (but not 
too low T), where dch/dfi ^ Ch/T and Xd oc Ch is ex- 
pected, and consequently axy oc —Bxd/T. Indeed, results 
for Nh = I indicate [|lj] that both crossovers Ts and T* ap- 
pear also as a change of sign in Xd{T) nearly at the same 
values. Here, the intermediate regime T* < T < Tg corre- 
sponds to an anomalous paramagnetic response Xd > 0. On 
the other hand, it is quite evident from our results that the re- 
lation is not valid at high T ^ t. Namely, in this regime 
al^ cx 1/T and a°y oc B/T'^ [0 is obtained, leading to 
R^{T ^oo) ^ const. On the other hand, from the high-T 
expansion Xd cx l/T"^ is acquired |]14|], so that the assumed 
relation would demand a^y oc B /T , in conflict with previ- 
ous a% (X B/T\ 

In conclusion, we have presented results for both dynamical 
and d.c. Hall constant within the t-J model on ladders and on 
square lattices. The main novel point is the observation of two 
crossover temperatures Tg and T* which are at low doping 
generally present in all systems. Both and are positive 
at T > Ts and change sign at Tg. While R^{T<Tg) stays 
negative, i?^ reveals a sign change into a hole-like behavior 
at r = r* < Tg as well as steep valuation of R^{T <T*). 
This reconciles some seemingly controversial theoretical re- 
sults [ pr] , |2|. Our results are in agreement with high-T ex- 
pansion results for R^{T) which at low c/j also show decreas- 
ing positive values with decreasing T. Quantum Monte Carlo 
results within the Hubbard model for Ru{iuj) correspond ef- 
fectively to high (imaginary) frequencies and low T, and be- 
ing negative they are in agreement with our findings for R^. 

How should we understand the above numerical results? 
At high T t and low doping Ch ^ 1, i?^ ^s well as i?^ 
are governed by a loop motion (that is where the dependence 
on i? 7^ comes from) of a hole within a single plaquette 
[jl^, [n|]. One expects R'^ > 0, but = 1/4 is a non- 
universal value which e.g. depends on the lattice coordination 
JTl]]. The electron-like r^{T^O) < represents an instanta- 
neous Hall response within the ground state near half filling is 



harder to explain, but is clearly the signature of strong corre- 
lations. On the other hand at low T, R^ tests the (low energy) 
quasiparticle properties. Evidently, at low doping and T < T* 
at least a single hole in an antiferromagnetic spin background 
behaves as a well defined hole-like quasiparticle leading to 

(T^O) - 1 both in 2D [|l7ll and in ladders [|]. Our re- 
sults for vanishing R^{T*) indicate that the quasiparticle 
character is essentially lost at quite low T ^ T* < J, with the 
pseudogap scale T*{ch) decreasing with doping. Such phe- 
nomenon is possibly consistent with the scenario of electrons 
being effectively composite particles (spinons and holons) in 
strongly correlated systems |Q, |T^, at least at T > T*{ch), 
whereby T*{ch) vanishes at optimum doping. 

Finally, let us note that our results for R^ are in several as- 
pects consistent with experiments on cuprates, and with LSCO 
in particular. At low doping c/i < 0.1 we find T* ~ J/2, close 
to the observed T* ~ 800 K. At the same time, we find a very 
steep dependence in R'^{T<T*). With increasing Ch, T*{ch) 
seems to have desired decreasing tendency, although to estab- 
lish this beyond a reasonable doubt more work is needed. 
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